arXiv:hep-ph/0112347vl 28 Dec 2001 


Some methods for the evaluation of complicated 

Feynman integrals 

A.V.Kotikov 

Bogoliubov Theoretical Physics Laboratory, Joint Institute for Nuclear Research 

141980 Dubna, Russia. 

Abstract 

We discuss a progress in calculations of Feynman integrals based on the Gegen- 
bauer Polynomial Technique and the Differential Equation Method. We demon¬ 
strate the results for a class of two-point two-loop diagrams and the evaluation of 
most complicated part of 0(l/iV 3 ) contributions to critical exponents of ^> 4 -theory. 
An illustration of the results obtained with help of above methods is considered. 


Last years there was an essential progress in calculations of Feynman integrals. It 
seems that most important results have been obtained for two-loop four-point massless 
Feynman diagrams: in on-shall case (see B 0) and for a class of off-shall legs (see ||). 
A review of the results can be found in [|j. Moreover, very recently results for a class of 
these diagrams have been obtained in the case when some propagators have a nonzero 
mass. 


In the paper, I review two methods for calculations of Feynman diagrams. 

The first one, so-called the Gegenbauer Polynomial Method (see [||] and also B3- 
||), has been used in particular for the evaluation of a s - corrections to the longitudinal 
structure function of deep inelastic scattering process. The structure of the results in 
Mcllin moment space (see 0-0) is very similar to the coefficients in [|I]. [J of the 
Mellin-Barnes transforms for the above doublc-bokses. The coefficients are similar also 
to ones which have arised (see [|14|]-[|T6j]) in expansions over the inversed mass for some 
two-loop two-point and three-point diagrams. 

A version of the second method, which is called as the Differential Equation Method 
I7|]-[|2(|, has been used in above calculations (see ||, £ij] and references therein). 

An illustration of some results which have been obtained with help of these two meth¬ 
ods is considered. The additional information about a modern progress in calculations of 


Feynman integrals can be found, for example, also in recent articles |2T], [22 


A. The Gegenbauer Polynomial Technique 


1 Basic Formulae 


Fifteen years ago the method based on the expansion of propagators in Gegenbauer series 
0) has been introduced in ®0. One has shown [|6j. || that by this method the 


see 


analytical evaluation of counterterms in the minimal subtraction scheme at the 4-loop 


1 





































level in any model and for any composite operator was indeed possible. The Gegenbauer 
Polynomial (GP) technique has been applied successfully for propagator-type Feynman 
diagrams (FD) in many calculations (see @,01). In the Section A we consider a develop¬ 
ment of the GP technique (obtained in |J), an illustration of the results obtained in 
and the application of the results calculated in (TT| . 


Throughout the Section A we use the following notation. The use of dimensional 
regularization is assumed. All the calculations are performed in the space of dimension 
D = 4 — 2e. Note that contrary to [§j we analyze FD directly in momentum ;r-space 
which allows us to avoid the appearance of Bessel functions. Because we consider here 
only propagator-type massless FD, we know their dependence on a single external mo¬ 
mentum beforehand. The point of interest is the coefficient function Cf, which depends 
on D — 4 — 2e and is a Laurent series in e. 

1.1 First of all, we present useful formulae to use of Gegenbauer polynomials. Fol¬ 
lowing 10- D-space integration can be represented in the form 


d D x = -Sr,_i(x 2 ) x dx 2 dx (A = D / 2 — 1), 


where x = f/Vr, and So- i = 27 T A+ 1 /r(A + 1) is the surface of the unit hypersphere in 
Rd • The Gegenbauer polynomials C^(t) are defined as |23], j7|j 


(1 - 2 rt + r 2 )-‘ = ]T C s n (t) r” (r < 1), Cj( 1) = 


n =0 


T(n + 25) 
n!T(2h) 


whence the expansion for the propagator is: 


(xi - x 2 ) 


25 


E 

n =0 


,XiX 2 ) 


l[xl 


2\ n / 2 


±L^e( x 2 2 -xl) + [xl^xfj 


where 


e(y) = 


1, if y > 0 


( 1 ) 


( 2 ) 


0, if y < 0 

Orthogonality of the Gegenbauer polynomials Cf(x) is expressed by the equation (see 


C A (fif 2 ) C^(x 2 x 3 )dx 2 = 


A 


c Cfo\x a), 


n + A 


( 3 ) 


where 8™ is the Kronecker symbol. 

The following formulae are useful (see BE3): 


c‘Jx) = E 

p>0 


(2x) n -' 2p (-l) p T(n - p + 8) 


(n — 2p)\p\T(8) 


(2x) 


n\ 




n—2p 


X, 


p> 0 


(n -2 p + 8)T(8) 
p\Y{n — p+8+1) 


and 


( 4 ) 
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Substituting the latter equation from ([|) for 5 = A to the first one, we have the 
following equation after the separate analyses at odd and even n: 



l"/2] / 

E C n-2pW- 


2 k + A)T(A) r (n + 5 - jfc)r(lfe + 5 - A) 
Jfc!r(<5) r(n + A + l-fc)r(<J-A) 


(5) 


1.2 Following || [nj we introduce the traceless product (TP) connected with 

the usual product x M1 ...x Mn by the following equations 


...fin _ 


X‘ 




g n!( l) p r(n p + A) gp 2p -iP 2 P x 2 p x p2 P +i x p n 


p> o 


*E 

p> 0 


2 2p p\(n — 2p)\Y(n + A) 

n!r(n — 2p + A + 1) 

(2 ) 2 Pp!(n — 2p)!r(n — p + A + 1) 


1/42 glJ-2p-lP2p j.2p X ^2p+1—Hn (g) 


Comparing Eqs.(ffl) and (13), we obtain the following relations between TP and GP 


^,fll.../l n ^,/Xl .../Z n _ 




n!r(A) c*(ii) (x^r' 2 , 


2 n Y(n + A) 


r(n + 2A)T(A) 2n 


2 r T(2A)T(n + A) 

We give also the simple but quite useful conditions: 


(7) 


...fin —— Z^ n == 


( 8 ) 


which follow immediately from the TP definition: x w-M»-w-Mn — q. 


The use of the TP xmakes it possible to ignore terms of the type that 
arise upon integration: they can be readily recovered from the general structure of the 
TP. Therefore, in the process of integration it is only necessary to follow the coefficient 
of the leading term x^ 1 ...x^ n . The rules to integrate FD containing TP can be found, for 
example in ]10], |Tl], [T^j. For a loop we have (hereafter Dx = (d D x)/(2 tt) d ) Q: 


Dx 


x■ 


Pi • "Pn 


All...fin 


x 2a (x 


y) 


2/3 


(4vr) D / 2 y 2(«+/3-a-i) 


nm A n ' a (u,P), 


(9) 


where 


A n,m (a, (3) = 


(ft )flm (/?) 


^?2+m(^ T fd A 1) 


and a n (a) = 


Y{D/2 — a + n) 

rR 


1 The Eq.(|)j) has been used in 00. 00 for calculations of the moments of structure functions of 

deep inelastic scattering. 
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Note that in our analysis it is necessary to consider more complicate cases of integra¬ 
tion, when the integrand contains 0 functions. Indeed, using the Eqs.(^) and (|7j), we can 
represent the propagator (x\ — x 2 )~ 2X into the following form 2 : 


1 _ 2’T(n + A) 

(xi — x 2 ) 2X “ rdr(A) 

Using the GP properties from previous subsection and the connection (|7|) between GP 
and TP, we obtain the rules for calculating FD with the 0-terms and TP. 


2(A+n) 


0(x 


X' 


! ) + 


_x 


(r 2 <- 

r 2 ^ 


X 2J 


( 10 ) 


1.3 The rules have the following form: 

y»/^l • • •f-l'n 


Dx 


x h 


x 2a (x — y) 2/3 


®{x 2 -y 2 ) 


) ,f-L 1 • • •fJ'n 


Dx 


rpf-Ll • ••fj'n 

d, _ , O 2> 


x 2a (x — y) 2 P 


E 


B(m, n\{3, A) 


(47t) d / 2 y 2 ( a+ 0 A P ttt. H-cxH-/5 — A — 1 


(/3=A) 




Q(y 2 -x 2 ) = 


( 47r )D/ 2 y 2 (a-i) r(A)(a-l)(n + A) 


( 11 ) 


1 y m-nn ^ B(m, n\/3, A) 


(47r) D / 2 ^(a+ZJ-A-i) ^ m 


09=A) 


T n — a + A + 1 
1 


(47r) D / 2 y 2 !"- 1 ) T(A) (n + A + 1 — a){n + A) 


( 12 ) 


where 


B(m, n\P, A) 


r (m + /3 + n) r(m + /3 — X) 

m!r(m + n + l +A)r(/3) T(/3 - A) 


The sum of above diagrams does not contain 0-terms and should reproduce Eq.(J|). 
To compare the r.li.s. of Eqs.(|ll],|T2|) and the r.li.s. of Eq.(9[) we use the transformation 
of 3 F 2 -hyper geo met r ic function with unit argument 3 F 2 (a, 6, c; e, b + 1; 1) (see |p5|): 


E 

fc =0 


r(fc + a)T(k + c) 1 _ r(a)r(l-a)r(6)r(c-6) 
k\T(k + f) k + b - r(/ — 6)r(l + b — a) 


r(i — a)T(a) 


E 


T(k + c - f + l)T(k + c) 


T(f - c)r(l + c - /) ^ k\T(k + 1 + c - a) k + c-b 


This is the case (when k = m, b = a + (3 — A — l,c = n + (3) to compare Eq.( 
sum of Eqs. ( 0 , 0 ). 

Analogously to Eqs.flTTD and (|T^ ) we have more complicate cases: 


(13) 
and the 


Dx 


x 


Ail-.-Mn 


x 2a (x — y) 


2/3 


0(x 2 


z = 


(47 r) D T 


y 


Atl-.-Mn 


e(y 2 


z 


-T-LA n \a,(3) 


y 


2(a+/3—A—1) 


2 In the case of the propagator (aq — x 2 ) 2S with <5 / A we should use also Eq. (§• 
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E 

m =0 


B(m,n\j3,\) / 

9 m 

Y y \ 

to 

1 

JO 

0 m+8+n 2 2\ 1 

fz 2 \ @(y 2 -z 2 ) \ 

g2(pt.-\-(3— A—1) \ 


m + a + /3 — A — 1 

\y 2 J m — a + n + l + XJ 


09=A) 


(4vr) D / 2 r(A) 


y 


fl\.. .[In 


_ Q(y 2 -z 2 ) 

y2(a—l) (q; — l){n + A + 1 — a) 


+ 


1 1 (Q(z 2 -y 2 ) fz 2 \ n+x Q(y 2 -Z 2 ) 


Dx 

OO 

E 

m=0 


X 


■ Ml • ■•fl'n 


- 2 (a-l) n _)_ \ 


<?>(z 2 -y 2 ) = 


x 2a (x — y) 2/3 
B(m,n\{3, A) / /y 2 


a — 1 

1 

(47r) D / 2 
Q(z 2 -y 2 ) 


y 


Ml • ••[^n 


n + 1 + A — a 
Q(z 2 -y 2 ) o 

y2{a+/3- A—1) ^ ( a >/5) 

' 2 ' m+/3+n 0(y 2 -z 2 ) 


(14) 


z 2(a+p- a— i) V\^ 2 / m + a + /3 —A —1 


C9=A) 


(4vr ) D / 2 r(A) 
1 1 
z 2 ( a_1 ) n + A 


[11... [l n 


Q{z 2 


\y A J m 

Q(z 2 -y 2 ) 


OL Tl -h 1 -f- A 


j/ 2 («-i) (q/ — l)(n + A + 1 — a) 


y 2 ) 


a 


2 . n +^ p 




0 ( 2 / 


,2 ^ 


n + 1 + A — a 


(15) 


One can easily see that the sum of the above diagrams lead to results identical to 


2 Calculation of complicated FD 


The aim of this section is to demonstrate the result of || for a class of master two-loop 
diagrams containing the vertex with two propagators having index 1 or A. 

Consider the following general diagram 


DxDy 


y 2a [z — y) 2t (z — x) 2 Px 2 z(x — y) 2s 


J(a,t,P, 7,s) 


and restrict ourselves to the FD A(a, (3, 7 ) = J(a, A, (3, 7 , A), which is the one of FD of 
interest for us here. It is easily shown (see [|26], [L 2 ], p|) that (a = 3 + A — (a + (3 + 7 )) 


C f [A(a,/3, 7 )] = C f [A(a,a,~f)\ = C f [J( 7 ,A,A,a,a)] = C/[J(a, 7 , A, A,/3)], (16) 

Doing Fourier transformation of both: the diagram A(a,/3, 7 ) and its solution in the form 
Cf[A(a, (3, 7 )](^ 2 ) _<t , where hereafter t — A + 1 — t, t = {a, (3, 7 , cr,...}, and considering 
the new diagram as one in the momentum x-space we obtain the relation 


C f [A(a,/3, 7 )] = 


ao(A)ao(a)ao(/3)a 0 (7) 

a 0 (S) 


C f [J(a, 1,^,7,!)] 


= ao(A)a 0 (a)a 0 (/l)ao(7)ao(c r ) C f [J(j3, 1, a, 7 ,1)] 


(17) 
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between the diagram, which contains the vertex with two propagators having the index 
A, and the similar diagram containing the vertex with two propagators having the index 
1 . 


Repeating the manipulations of j |2fi| . [27]. [12|. [J we can obtain the following relations: 


C f [J((3,l,a,7,l)\ = C f [J0, l,cr, 7 ,1)] = Cf[J(l,l, 7 ,a,a)] = C f [J(a, 1,1, 7 , /3)] (18) 

Thus, we have obtained the relations between all diagrams from the class introduced in 
the beginning of this section. Hence, it is necessary to find the solution for one of them. 
We prefer to analyze the diagram A(a,/3, 7 ), that is the content of the next subsection. 

2.1 We calculate the diagram A(a,/3, 7 ) by the following way * 3 4 * * : 


. . ( 10 ) 2 7 T(n + A) 

A{a '^ = E^p 


DxDy 


yf-ll •••H-n 




71 =0 


x 2 ^(z — x) 213 y 2a {x — y ) 


2A 


G(z 2 -y 2 ) + 0{y 2 -z 2 ) 


^2(n+A) 


y2(n-\-X) 


P4M) 


(4tt) d /z r(A) a - 1 ^ ri!r(A) 


E 


2’T (n + A) 


Dx 


gl-Ll • ••fJ'Tl ••l-l'Tl 

x 2 ^(z — x) 2 h 


G(z 2 - x 2 ) 


0(x 2 - z 2 ) 


A + n + 1 — a y z 2 ( n+ Ax 2 ( a b x 2(n+\) z 2(a 1 ) 


_ 1 Q(z 2 -x 2 ) Q(x 2 — z 2 ) ' 

A + n + a — 1 ~ 2 (n+A+a— 1 ) ^(n+A+a-l) 

After some algebra we have got (see ||) the result in the form: 

111 


C f [A{a,M)] = 


(An) D T(A) a - 1 L 


(19) 


where [] 


/ 


T (n + 2A) 

1 

A + n + a — 1 


—————-- (A n,0 (a — 1 + 7 , (3) + A n ’°(n + A + 7 , /3)) 

A + n + 1 — a \ / 

(aT’ 0 ( 7 , P) + A n -°(n + 0 + A + 7 - 1 ,P)) 


( 20 ) 


(n) 

3 The symbol = marks the fact that the equation (n) is used on this step. 

4 We would like to note that the coefficients in Eqs.fg^) and (21) are similar to ones (see [^5|) appeared 

in calculations of FD with massive propagators having the mass to. The representation of the results for 
these diagrams in the form ‘finiz 2 /m 2 ) n (ip n are the coefficients, which are similar to ones in Eqs.([20|) 

and (ml)) is very convenient to obtain the results for more complicated FD by integration in respect of 


to (see ]17|- [20f) of results less complicated FD. 
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~ r(l - (3)T(\ + 1 - a)r(A - 1 + a)r(l ~/3 + A)r(l - 7 )r(a + /3 + 7 - A- 2 ) 
r( 2 A)r (2 + A - a - f3)T(a + 7 - 1 )T (2 + A - 7 - f3)T(a + (3 - A - 1) 

V r(n + 2A) HgtM 1 f21 x 

^ n!T( 2 A) r(^-A) A + n + a -1 1 J 

r(a + (3 + 7 — A — 2 )T (2 — a — 7 ) 

T(3 — a — /3 — 7 — n)T(a + 7 + A — 1 + n) 

r(i- 7 )rQ3 + 7 -A-i) 

+ r(7-A-n)r(2 + 2A-/3-7 + n). 

Thus, a quite simple solution for A(a,/3, 7 ) is obtainedf]. In next section we will con¬ 
sider the important special case of these results. 

2.2 As a simple but important example to apply these results we consider the diagram 
J(l, 1,1,1, a). It arises in the framework of a number of calculations (see |2!| [3U], [H], 
|32|). Its coefficient function 1(a) = C/[J(1,1,1, l,a)] can be found (see |J) as follows 


1(a) 


g/WA.A,A,A, a) ] and 
Cf[J( A, A, A, A, a)) = Cf[A(a, 3 — A — a, A)] 


From Eqs. (|2D|) and ([2T|) we obtain 


1(a) 


X 


2 r 2 (A)r(A-a)r(a + l-2A) 

~ (in) D r(2A)r(3A — a — 1) 
T 2 (l/2)r(3A - a - l)r( 2 A - a)r(a + 1 - 2 A) 
_ r(A)r(2A + 1/2 - a)r(l/2 - 2A + a) 


OO 

+ E 


n =0 


T(n + 2A) 1 

r(n + a + l) n + 1 — A + a 


Note that in Kazakov has got another result for 1(a) 


( 22 ) 


2 r 2 (A)r(l-A)r(A-a)r(a + l-2A) 
Mi (4n) D r(2A)r(a)r(3A-a - 1) 

T(A)r (2 - A)r(a)r(3A -a- 1) 
r(2A- l)r(3-2A) 

_ y^/ Nn r(n + 2A) /_1_ 

J r(n + 2 — A) Vn + l-A + a 

—fl x ' 


n + 2A 


a 


(23) 


5 Before our studies, the possibility to represent Cf [A(a, /3, 7 )] as a combination of 3 -F 2 -hypergeometric 
functions with unit argument, has been observed in psf. 
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From Eqs. (p2|) and (|23|) we obtain the transformation rule for 3 F 2 -hypergeometric 
function with argument — 1 : 


3 .F 2 (2a, b, 1; b + 1 , 2 — a; —1) 


T(2 — a)r(6 + a — 1)T(& — a)r(l + a 
r(2a)r(l + 6 — 2a) 


7 —;- 7 ' 3-^2 ( 2 a, b, 1 ; b + 1 , 6 + a; 1 ) 

o + a — 1 

-— -- ■ 3 F 2 (2a , 1 + a — 6 ,1; 2 + a — 6, 2 — a; — 1), 

1 + a — b 



where a = A and b = 1 — A + a are used. 

Equation (^4|) has been explicitly checked at a = 1 and b = 2 — a (i.e. A = 1 and 
a = 1), where the 3 F 2 -hypergeometric functions may be calculated exactly. It is very 
difficult to prove Eq.(^4|) at arbitrary a and b values: the general proof seems to be non¬ 
trivial. Note that it is different from the equations of [33] and may be considered as 
a new transformation rule. 


3 Applications 


3.1 The above results have been used for evaluation of very complicated FD which 
contribute mostly in calculations based on various type of 1/IV expansions: 


In the calculation (in |54]) of the next-to-leading (NLO) corrections to the value of 
dynamical mass generation (see |33]]) in the framework of three-dimensional Quan¬ 
tum Electrodynamics. 


• In the evaluation (in |3f|) of the correct value of of the leading order contribution 
to the /3-function of the 0-term in Chern-Simons theory. The /3-function is zero 
in the framework of usual perturbation theory but it takes nonzero values in 1/N 
expansion (see [ |37jj ). 

• In the evaluation (in |32|]) of NLO corrections to the value of gluon Regge trajectory 
(see discussions in j32] and references therein). 


In the calculation (in |3£|) of the next-to-leading corrections to the BFKL intercept 


of spin-dependent part of high-energy asymptotics of hadron-hadron cross-sections. 


In the calculation (in |38], [hj) of the next-to-leading corrections to the BFKL equa¬ 
tion at arbitrary conformal spin. 

In the evaluation (in [Q) of the most complicated parts of 0(l/iV 3 ) contributions 
to critical exponents of </> 4 -theory, for any spacetime dimensionality D. 


We consider here only basic steps of the last analysis ]2^| . Since the pioneering work 
of the St Petersburg group |26|. [hj, exploiting conformal invariance [^l| of critical phe¬ 
nomena, it was known that the 0(l/iV 3 ) terms in the large-iV critical exponents of the 
























non-linear cr-model, or equivalently 0 4 -theory, in any number D of spacetime dimensions, 
derives its maximal complexity from a single Feynman integral /(A) (see [[ff]]): 


/(A) 


d 

dA 


lnII(A, A) 


A=0 


(25) 


where 


~>2(A+A) 

n(A,A) = — 7 — 


d D yd D z 

y 2 z 2 [x — y) 2X (x — z) 2X (y — z) 2 ( A+A ) 


(26) 


is a two-loop two-point integral, with three dressed propagators, made dimensionless by 
the appropriate power of x 2 . 

The result, obtained by GP technique, is 


/(A) = T(l) - T(1 - A) + 
where d>(x) = r'(x)/r(x) [] and 


$(A) - §tt"(A) - £*"(1) 
T'(l) - T'(A) 


(27) 


2A-1 


4(A) =4/ dx ^— — j l.ijj —- l.iA-1 )} (Li 2 (a:) = J^V/n 2 ), 


(28) 


n> 0 


In [ [43| , |24| the integral /(A) has been expanded near D = 2 and D = 3, respectively, 
(i.e. for D = 2 — 2e and D = 3 — 2e) up to e 8 in the form of alternative and non-alternative 
double Euler sums [011 0H1. 


3.2 The Eq. 


together with the “uniqueness” relations [[26) , [27| , |29j and the 


integration by parts extended in [|1CJ-||T2[ for the form of massless propagators 

~ £/u,. -./An /( x 2 y i ias been used for the evaluation of the aq-cor reef ions to the longitudinal 
structure function Fl of deep-inelastic scattering process [|. The corresponding results 
1 contain the sums 


10|, 113 


/ \ 1 n (—1 i m+1 

*>(«) = (1 - (-ir)-« 2 ) + (-ir 57 ^-^ 5 —. 

m= 1 

/ N 3 n (— 1 i m + 1 

*a(n) = (l-(,!)")-((SJ + f-ir^i—L-, 

m= 1 

C7 ^ / 1 \ T7T-—|— 1 

A'wW = (l-(-l)”) 5 C(3) + (-l)”J]^5— 5i(m) (29) 

771=1 

6 The function $(A) is quite similar to most complicated part of the NLO corrections to BFKL intercept 
(see |f2|, [is|, |39|] and references therein). 

' In the framework of supersymmertic extension the evaluation of gluino contribution to Fl has been 
done in Q . 

8 We would like to note that the results contain an error in gluon sector, which is not essential 

for phenomenology. The correct results have been found in 
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which can be obtained by direct calculations (with help of the optical theorem 9 ) only at 
even n. The results for the anomalous dimensions of Wilson operators contain also the 
function (p9|) (see |5T[| ) . 

The analytical continuation of the results (|29|) to integer values and to real ones (and 
even to complex ones) can be found, respectively, in |11[ and . The continuation to the 


integer values has been wide used in fits of deep inelastic experimental data at the next- 
to-leading-order (NLO) approximation [53, 54, |55[ and at the next-next-to-leading-order 
(NNLO) level (5(|-|)8|]. The continuation to the real values is very important for small 


Bjorken x phenomenology. In the Ref. |52[] the extension of previous results |60], |6T| has 
been performed for an approximation of Mcllin convolution by a sum of usual products. 
The extension give a possibility to obtain the following results: 


• To extend (in |62|) the solution of DGLAP equation in double-logarithmic approx¬ 
imation (see |63|) to NLO approximation. 

• To explanate (in [|64[]) a sharp change of the “intercept” value ap at Q 2 ~ 2 GeV 2 
for the power-like asymptotics of parton distributions and structure function F 2 , 
observed in [p5fl. 


To demonstrate (in 0) the positivity of the value for longitudinal structure func¬ 


tion of deep-inelastic scattering at small x range. The results have been obtained in 
so-called renormalization-invariant perturbation theory (see Ref. [[54| and references 


therein), which resums properly the large and negative NLO corrections (see |67 
in the gluon part of the longitudinal Wilson coefficient function. 


• To extract at small x values the gluon density in |68|] (following to the articles ||69l] ) 
and the longitudinal structure function F Tj (in ]7D[) from the experimental data for 
the structure function F 2 and its derivation dF 2 /d\nQ 2 . 


B. The Differential Equation Method. 


The idea of the Differential Equation Method (DEM) (see [|L7[-[|T9| and reviews in 
|): to apply the integration by parts procedure |7T], |6]] to an internal n-point subgraph 
of a complicated Feynman diagram and later to represent new complicated diagrams, 
obtained here, as derivatives in respect of corresponding masses of the initial diagram. 

The integration by parts procedure (see also |L7[]-[^0|) for a general n-point 

(sub)graph with masses of its lines mi,m 2 , ...,m n , line momenta pi,p 2 — p\ — pi 2 ,p n = 
Pi ~ Pin and indices R, R, - -, jR, respectively, has the following form: 


0 = 


d D pi 


d 

dpt 


-t 


A 4 

Pi 


LD 


i =1 


(30) 


9 The optical theorem is very powerful in calculations with massive particles, too (see j50| and references 
therein). 
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-1 


d °pi (n ^ 


i=l 


D - 2 3l 


1 _ U!± 

Cl 


n 


i=2 


ji 


m\ + ml + pfj — ci 


where Ck = pi + ml are the propagators of n-point (sub)graph. 

Because the diagram with the index (ji + 1) of the propagator q may be represented as 
the derivative (on the mass Wj), Eq. (j3Gp leads to the differential equations (in principle, to 
partial differential equations) for the initial diagram (having the index ji, respectively). 
This approach which is based on the Eq.(|30|) and allows to construct the (differential) 
relations between diagrams has been named as Differential Equations Method (DEM). For 
most interested cases (where the number of the masses is limited) these partial differential 
equations may be represented through original differential equation^, which is usually 
simpler to analyze. 

Thus, we have got the differential equations for the initial diagram. The inhomo¬ 
geneous term contains only more simpler diagrams. These simpler diagrams have more 


trivial topological structure and/or less number of loops [[L7|] and/or ends [[K 

Applying the procedure several times, we will able to represent complicated Feynman 
integrals and their derivatives (in respect of internal masses) through a set of quite simple 
well-known diagrams. Then, the results for the complicated FD can be obtained by 
integration several times of the known results for corresponding simple diagrams Q. 

Sometimes it is useful (see [[FfJ) to use external momenta (or some their functions) 
but not masses as parameters of integration. 


The recent progress in calculation of Feynman integrals with help of the 
DEM. 


1. The articles 14 and 15 


a) The set of two-point two-loop FD with one- and two-mass thresholds has been 
evaluated by DEM (see Fig.l). The results are given on pages 2 and 3 and of some of 
them have been known before (see ||T4[). The check of the results has been done by Veretin 
programs (see discussions in [14]. |l(J and references therein). 


b) The set of three-point two-loop FD with one- and two-mass thresholds has been 
evaluated (the results of some of them has been known before (see |T1|])) by a combination 
of DEM and Veretin programs for calculation of first terms of FD small-moment expan¬ 
sion (see discussions in ||T4|, p~H|| and references therein). 


2. The article 75 : 


10 The example of the direct application of the partial differential equation may be found in 172 . 

11 In calculations of real processes (essentially in the framework of Standard Model) it is useful to use 
the relation (1) (at least, at first steps of calculations) to decrease the number of contributed diagrams 
(see |l^|-[|T9) and [73j and references therein). 
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Figure 1: Two-loop self-energy diagrams. Solid lines denote propagators with the mass 
m; dashed lines denote massless propagators. 


The full set of two-point two-loop on-shell master diagrams has been evaluated by 
DEM. The check of the results has been done by Kalmykov programs (see discussions in 
0 |76| and references therein). 


3. The articles 


The set of three-point and four-point two-loop massless FD has been evaluated. 


Here we demonstrate the results of FD are displayed on Fig.l. 

We introduce the notation for polylogarithmic functions |77| : 

t; _ c ^ c f„\ _ (-!) a+b f 1 log a (t)log b (l -zt) ^ 

Li a\Z) S a — i i yZj , S a -\-i byZ) , , , / dt. 

a! b\ J 0 t 

We introduce also the following two variables 

f 

771 


q- i - V z /( z ~ 4 ) 
z = —, y = 


1 + y/z/iz ~ 4) 


ThenQ 


12 We would like to note that the coefficients of expansions of the results (0) in respect of z are very 
similar (se e |lo| -II) to results for the moments of structure functions of deep inelastic scattering, i.e. to 
the sums (] 2 ^) . 
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q 2 ■ h = ~ log 2 (-^) log(l - z) - 2 log(-z)Li 2 (z) + 3Li 3 (z) - 6 Si, 2 (z) 

-log(l-z)(c 2 + 2 Li 2 (z)) , 

q 2 ■ h = 2(2 log(l + z) + 2 log(—z)Li 2 (—z) + log 2 (-^) log(l + z) + 4 log(l + z)IA 2 {z) 

— 2Li 3 (— z) — 2Li 3 (z) + 2 S'i i2 (^ 2 ) — AS\p{z) — 4S'i j 2 (— z), 

q 2 ■ 1 12 = Li 3 (z) — 6(^3 — C 2 l°g U lo S 3 y “ 4 log y Li 2 (y) 

o 

+ 4Li 3 (y) - 3Li 3 (-y) + ^Li 3 (-y 3 ), 
q 2 ■ 113 = —6Si i2 (z) - 2 log(l - z) (^2 + Li 2 (z)) , 

q 2 ■ 1 14 = log(2 - z) ^log 2 (1 - z) -2 log(-z) log(l - z) - 2Li 2 (z)^ 

~ | log 3 (l - z) - 2 C 2 log(l - z) 

+ log(-z) log 2 (l-^) 

— * 51 , 2 ( 1 /(1 — ^) 2 )+2*5i i 2 (1/(1 — z))+2Si t2 ( — 1/(1 — z))+- log 3 y 
+ log 2 y (2 log(l + y 2 ) - 3 log(l -y + y 2 )) 

- 6(3 - Li 3 (— 2 / 2 ) + ^Li 3 (— 2 / 3 ) - 6 Li 3 (—y) 

+ 2logy ^Li 2 (—y 2 ) - Li 2 (—y 3 ) + 3Li 2 (-y)) , 

q 2 ■ hb = 2Li 3 (z) - log(-z) Li 2 {z) + C 2 log(l - z) 

+ \ log 2 y (8 log(l - y) - 3 log(l - y + y 2 )) - 6(3 

+ l log 3 y - ^Li 3 (—y 3 ) + 3Li 3 (—y) + 8 Li 3 (y) 
o 3 

+ logy (Li 2 (—y 3 ) - 3Li 2 (—y) - 8 Li 2 (y)j , 

q 2 • I 123 = -C 2 (log(l -z)+ logy) - 6(3 - ^ log(l - y + y 2 ) log 2 y + Li 3 (-y 3 ) - 9Li 3 (-y) 

- 2 logy ^Li 2 (—y 3 ) - 3Li 2 (-y) j , 

q 2 ■ I 125 = -2 log 2 y log(l-y) - 6(3 + 6 Li 3 (y) - 6 logy Li 2 (y), 

? 2 ■ 5 1234 = — 6/3 — 12Li 3 (y) — 24Li 3 (—y) 
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+ 8 log y (h\ 2 (y) + 2Li 2 (-2/)) +2 log 2 y (log(l - y) 
+ 2 log(l + y)^ . 


( 31 ) 


Here we demonstrate the results of FD are displayed on Fig.2. 

We consider here the following master-integrals in Euclidean space-time with dimen¬ 
sion D = 4 — 2e: 


ONS{lJ}(i,j,m) 

= A" 1 J d D kP^(k,lm) P U) (k - p, , 

, k , m) 

= A" 2 J i D k 1 d D knP^{k,,lm)P ii \k 1 -k 2 ,Jm) 


x P {k \k 2 -p,ICm )\ p2= _ m2 , 

V{lJJC£}{i,j,k,l,m) 

= K~ 2 J d D kid D k 2 P^ {k 2 - p, 1m) 


xPWfa - k 2 ,Jm)P ( - k \k 1 ,Km) P^(k 2 ,Cm)\ p2 

~F{ABTJ1C\(a , b, i,j, k , m) 

= m 2 K~ 2 J d D k 1 d D k 2 P (a \k 1 ,Am)P w (k 2 ,Bm ) 


xP (i) (fci -p,lm) 

xPV\k 2 -p,Jm ) P^\k, - k 2 ,JCm) | 2 , 


where 


K = 


r(i + e) 


(4n 


2 [m 2 


P ( ' l \k , m) = 


(k 2 + m 2 ) 1 ’ 


the normalization factor 1/(27t) d for each loop is assumed, and A, 13,1, J J C = 0,1. 

The finite part of most of the F-type master-integrals can be obtained from results of 
Ref.[n] in the limit z 1. F10101 and Flllll have been calculated in Refs. |78|, |75[| , 
respectively. Instead of the usually taken F01101 integral [[78], |80]] we consider Jill as 
master integral. We recall the results of all master integrals for completeness. The last 
master integral F00111 has been found in |75|| . 

The finite part of the integrals of V- and J-type can be found in Refs.[^l|]. The 
calculation of some ~ e and parts of master integrals of this type have been performed 

by DEM. 

The results for F-type master-integrals are follows: 


7T 


¥{ABIJK}{\, 1,1,1,1, m) = ai C(3) + a„S 2 + a z inC,(2) + 0{e), (32) 

v 3 

and the coefficients {a*} are given in Table I: 
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TABLE I 


Flllll 

F00111 

F10101 

F10110 

F01100 

F00101 

F10100 

F00001 

«i 

1 

0 

-4 

-1 

0 

-3 

-2 

-3 

a 2 

y 

2 

9 

27 

2 

9 

27 

2 

27 

2 

9 

0 

a 3 

0 

0 

1 

3 

0 

1 

1 

2 

3 

1 


where |82|, |78|. 1771 


S 2 = -^=C1 2 = 0.260434137632 ■ ■ ■. 

9VS V3/ 

Here we used the m 2 — ie prescription. The results for the remaining master integrals are 
the following ones: 


VdJ^Kl. 1 , 1 , l.rn) = T + i (£ _ JL) + h + | c(2) _ 4^ - “ft 


7T 

7 ! 


In 3 + £ 


Hi 7 T Q nr 

— + b 2 ((2) — b 3 ((3) — 12—^= — 635^ + & 4 C( 2 ) In 3 + - 64 ^—^= 


63^ 7 r 1 7T 9 65 7T 2 lLs 3 (^f) 

—5, In 3 + 4— In 3- ■= In 2 3 - —-=C(2)- 

4 2 2^3 2 V ; 2 v/3 


+ 0(e 2 ), 


where the coefficients {fcj} are listen in Table 11 
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TABLE II 


b\ 

b 2 

b 3 

bi 

b§ 

Villi 

-1 

-6 

y 

2 

4 

9 

V 1001 

3 

8 

3 

2 

0 

21 


Jlll(l, 1,1, m) 


J011(l, 1, 2, m) 


17 59 




m 2 —- H-1-h e< -h 8C(2) 

V2e 2 4e 8 I 16 ; 

, f1117 


65 


i 


32 

1 -As 


52C(2) + 48C(2) In 2 - 28C(3) \ + 0(e 


2(1 — 2 £ )(1 - 3e) (e 2 + 2 V3 3^ 2 * 

'•j= - §C(2) - 6 -^= In 3 + |C(3) + 27 S 2 }• + 0(e ‘) |, 


(33) 


( 34 ) 


13 The results for the master integral V1001 had a little error (see |8^|) 
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J011(l, 1,1, m ) 


m 2 4 — 15e: / 1 3 n 

'~2 (1 - 2e)(l - 3e)(2 - 3e) V? + 27! 
. 45 7T 9 7T 81 

+ ^T^- 2 VS ln3 + T & 


+ £ 2 |l 2 -C( 2 )-^!s 2 - 

27 7T , 9 7T .. . 

- "i" vf 113 “ 21 vS ^ 2 ” 


867 7r 207 7r 


+ 


243 


32 v/3 

si mu 

2 y/3 


h- In 3 4-— S 2 In 3 

V3 4 


0{£ 3 


ONSll(l, 1, m) = 


1 - 2c 


7T 


V3 “ \y/3 


71 


- 7 = + £ < In 3 — 95*; 


(35) 

(36) 


+ .i' ) 9.S'i In 3 - ln -3 -6—lL4 

2 y/3 V3 


7 r 


- 3 ^jC( 2 ) > + °( £3 ) 


where 77 


/ In 2 

9 

2 sin - 

Jo 

2 


d6 and Ls, I — ] = -2.14476721256949 • • • 


Ls 3 (x) = 


The above results were checked numerically. Pade approximants were calculated from 
the small momentum Taylor expansion of the diagrams |84|] . The Taylor coefficients were 
obtained by means of the package |S3| with the master integrals taken from Further 
we made use of the idea of Broadhurst f87| | to apply the FORTRAN program PSLQ |[77 


to express the obtained numerical values in terms of a ‘basis’ of irrational numbers, which 
were predicted by DEM. 

Let us point out that the numbers we obtain are related to polylogarithms at the sixth 
root of unityQ and hence are in the same class of transcendentals obtained by Broadhurst 
|87H in his investigation of three-loop diagrams which correspond to a closure of the two- 
loop topologies considered here. 
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Figure 2: The full set of two-loop self-energies diagrams with one mass. Bold and thin 
lines correspond to the mass and massless propagators, respectively. 
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